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FIGURE la 

(b+d)+i(a+c) 




FIGURE lb 




(d-b)+i(c-a) 
(d-b)-i(c-a) 

FIGURE lc 
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FIGURE 2 
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Routine: main 

Data, Variables, and Effects: 

input: data presented as a function f(p,q), 
giving the function value at the 
point (p-iq)/(p+iq) 
output: array of complex Fourier coefficients 
effect: calculates Fourier transform 



1. Normalize input: 

a) Initialize arrray of complex fourier coefficients to zero 

b) Compute abar, bbar, cbar so that 

p^-cp" -2 pq 

f(p,q)=f(p.q> Cabar* — + bbar*— - — +cbar] 

p^+q 2 p 
is zero for (p,q)=(l,0), (0,1), and (1,1) 
Then: f(p,q) returns normalized values above 



2. Start recursive calculation of arrows on top of 
the circle with call to gener (Figure 6) withinitial 
values: arrow= doe, generation=l,envy=Q, f(-l,l) 



3. Start recursive calculation of arrows on bottom of 
the circle with call to gener (Figure 6) withinitial 
values: arrow= doe, generation=l,en vy=Q, f(l,l)=0 

,.t , 

|4. Normalize output: divide Fourier coefficients by pi } 



1 5. Display computed coefficients and exit 



FIGURE 5 
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15. End cascade and make 
correction calling prune 
(Figure 8) with arrow as 
argument 



envy >= NVAN and 
.generation >= MINGy 



± 



7. Generate two descendant 
arrows in cascade using 
least-squares regularization 



recursion 
x2 



8. Update Fourier coefficients 
calling charles (Figure 7) 
with each descendant arrow 



14. Recursively call gener 
with the two descendant 
arrows, envy <-- envy, 
generation<~generation+ 1 

Routine: gener 

Data, Variables, and Effects: 

input: arrow, generation, envy, function value 
output: void 
effect: recursive routine to generate arrows and 
update Fourier coefficients 




FIGURE 6 
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1 6. Clear flag I 







17. Calculate preliminary 
quantitites 



Routine: charles 

Data, Variables, and Effects: 

input: one arrow 

output: flag representing whether array 

update for arrow is non-negligible 

effect: modifies array of Fourier coefficients 

in accordance with the formulas in Step 2 
for the Fourier transform 



19. Pass data to routine for 0, 
+1,-1 Fourier coefficients 



27. Proceed to next 
Fourier coefficient 



25. 

Is contribution 
^non-negligible^/ 




21. Set flag 





N 




24. Compute contribution of 
current arrow to current 
Fourier coefficient 







FIGURE 7 
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28. Calculate final updat 
trigonometric functions 
the end of the discussion 


es taul,tau2 of the two lagged 
For terminal processing as at 
of Fourier transform 







for terminal processing as at the end of the discussion 
of Fourier transform using 3 further data samples for each of 
the two descendant arrows 



Routine: prune 

Data, Variables, and Effects: 

input: one arrow 
output: void 

effect: modifies array of Fourier coefficients 

by implementing the final processing of 
arrows in the cascade 



30. Modify the Fourier coefficients by adding those of 
sigmal-taul and sigma2-tau2 truncated as at the end 
of the discussion of Fourier transform 




FIGURE 8 
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34. Calculate modified Fourier coefficients 

n=£0,-l,+l n 
for m=0,-l,+l, where the coefficients c m 
are given at the beginning of the section 
on inverse Fourier transform 



35. Using the formulas from Step 1 at 
the beginning of the section on inverse 
Fourier transform, calculate the wavelet 

coefficients e A for A = I, U, T, U' 1 , T _1 , 

U ~ 2 T ' 2 U A T l , T 1 U' 1 



36. For A=U, T, T ' 2 U" 1 T" 1 , T _1 U "\ IT 2 , 
use the wavelet coefficients from program 
segment 2 to stuff an arrow labeled A 
with the correct lagged trigonometric 
function and wavelet coefficient and 
serially call gener (Figure 10) with these 
arrows as argument thus starting these 
six recursive calculations of arrows 



Routine: main 

Data, Variables, Effects: 

input: array of complex Fourier coefficients 
c for In! <= N 

output: array of three-tuples (p,q,z), where 
output function takes complex value 
z at (p-iq)/(p+iq) given in the counter- 
clockwise order on the circle 

effect: calculates inverse Fourier transform 



37. Correct the outp 
adding Cq + c^ e 


ut function values by 
-inG + <£j ein© 






38. Display computed function values 
and exit 



FIGURE 9 
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recursion 
x2 ' 



41. Recursively call gener with the 
two descendant arrows as argument 



39. Generate the two descendant arrows 
in the cascade calculating arithmetic 
wavelet coefficients using the formulas 
from Step 1 from the section on inverse 
Fourier transform 




Routine: gener 

Data, Variables, and Effects: 



input: one arrow 
output: void 

effect: recursive routine to generate 
arrows and compute function 
values 



42. Perform the final update of 
lagged trigonometric functions 



43. Stuff output array with function values 
at endpoints of the two descendant arrows 



FIGURE 10 



